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A given dynamics for a composite quantum system can exhibit several distinct properties for 
the asymptotic entanglement behavior, like entanglement sudden death, asymptotic death of entan- 
glement, sudden birth of entanglement, etc. A classification of the possible situations was given in 
[M. O. Terra Cunha, New J. Phys 9, 237 (2007)] but for some classes there were no known examples. 
In this work we give a better classification for the possibile relaxing dynamics at the light of the 
geometry of their set of asymptotic states and give explicit examples for all the classes. Although 
the classification is completely general, in the search of examples it is sufficient to use two qubits 
with dynamics given by differential equations in Lindblad form (some of them non-autonomous). 
We also investigate, in each case, the probabilities to find each possible behavior for random initial 
states. 
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I. INTRODUCTION 

Entanglement is a fundamental property of compos- 
ite quantum systems, first noted by Schrodinger p]]. The 
best knowledge of the whole of a composite quantum sys- 
tem may not include complete knowledge of its parts. It 
has strong conceptual implications on physics, since it is 
a property that has no classical analog, so we are forced 
to change significantly our perspective of Nature. Such 
peculiar character allows it to be considered as a funda- 
mental resource for some non-classical tasks as telepor- 
tation of a quantum state quantum computation 
quantum cryptography Q, etc [2^1 • Once entanglement 
is considered a resource it seems natural to quantify it 
In all the applications named above, it is necessary 
to optimize the amount of entanglement in a suitable 
composite quantum system to best execute the desired 
task. 

Real quantum systems always interact with its envi- 
ronment, irrespectively of the efforts to protect it. This 
interaction will, in general, create some entanglement be- 
tween the quantum system and the environment, and this 
entanglement will, somewhat ironically, spoil the entan- 
glement between the parts of the "useful" system (for 
bipartite systems, this affirmation has a precise meaning 
provided by the monogamy of entanglement theorem [3]). 

While in most of the models used to describe quantum 
open systems the coherences of a state decays asymptot- 
ically to zero, it was recently recognized that entangle- 
ment may "die" at finite time a phenomenon called 
entanglement sudden death This phenomenon has 
called some attention, specially connected to the diffi- 
culty of keeping entanglement alive for its uses as a re- 
source. Some interesting generalizations were studied 
[lo| . and some experiments were proposed [ll| and re- 
alized This phenomena, though, has a simple ex- 
planation if one looks at the geometry of quantum states 
Namely, while the set of "decohered" states always 



have zero volume inside the set of all possible quantum 
states, the set of separable states has not only a positive 
volume but also non-empty interior when the global 
system have a finite dimensional Hilbert space. 

The geometrical approach to the problem allows one 
to classify the dynamics of a quantum system according 
to the geometry of its asymptotic states (if the dynamics 
implies them) relative to the set of separable states 
In the cited paper some classes were exemplified, but to 
that time it was not clear whether all a priori possible 
situations could be found. 

In this paper we review the geometric classification of 
entanglement dynamics and provide explicit examples to 
all a priori possible situations. All examples are given 
in the two-qubit Lindblad differential equations context, 
with some cases using non-autonomous equations (ex- 
actly those in the classes for which examples were not 
previously known). We also introduce a new analysis of 
how often each specific behavior occur for a given dy- 
namics, in the light of probability theory applied to the 
set of initial states [Tsj . 



II. THE GEOMETRY OF ENTANGLEMENT 
SUDDEN DEATH: GENERAL PICTURE 

What can we say about the geometry of entanglement, 
or the geometry of the set of separable states, for general 
multipartite systems? First of all, that the set of sep- 
arable states is closed, convex and with non-empty in- 
terior (we shall assume finite dimensional Hilbert spaces 
throughout the paper). Its complement relative to the set 
of quantum states also has non-empty interior and is cer- 
tainly non-convex. Actually, in general, its complement 
is much larger, i.e., it has greater volume (if one consider 
the Hilbert-Schmidt metric, for instance). An extremely 
oversimplified illustration of this situation is given in Fig. 
[TJ We call here D the set of all quantum states and we 
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are going to consider it immersed in the set A of Hermi- 
tian matrices of unity trace; S the subset composed by 
the separable states, dS and dD their boundaries relative 
to D and A, respectively; E — D — S the set of entan- 
gled states. The boundary of the set of quantum states 
is composed by all that states which have at least one 
zero eigenvalue so, in particular, it contains all the pure 
states. Note that there are both entangled and separable 
pure states in dD. Actually, more than that, the "area" 
of the separable states inside dD is non-zero 




FIG. 1: Diagram of the set of entangled states. 

Let us consider a dynamics with a non-trivial station- 
ary set, St. By stationary set we mean that for every 
initial state p and open set V ^ St we have that p{t) 
(the state at time t) belongs to V for all t sufficiently 
large. Of course, if some dynamics accepts a set of sta- 
tionary states, this set will be the smallest stationary 
set of the dynamics. Anyway, from the simple picture 
given in Fig. [H and considering the location of St in it, 
we may distinguish three possibilities which have conse- 
quences to the asymptotic dynamics of entanglement: i) 
St C Int(S') implies that every initial entangled state 
will lose all of its entanglement at finite time (sudden 
death of entanglement) ; ii) \iStr\dS ^ 0, than, only with 
that information, many situations can occur: asymptotic 
or sudden death of entanglement and non-zero asymp- 
totic entanglement; iii) if St C E, every initial state ex- 
hibit some entanglement asymptotically. 

The complete classification must yet consider that the 
stationary set, St, can consist of a single state {e.g., ther- 
mal equilibrium state) or by a non-trivial set {e.g., for 
phase reservoirs). In this sense, each situation above 
gives rise to two cases, in a total of six classes. 

Note that, in cases ii) and iii), if we start with a sep- 
arable state it is possible in the first one and certain in 
the second, that entanglement will be created, a situation 
which may be called sudden birth of entanglement It 
is good to stress that, since the only information we have 
about the dynamics is some partial information about a 
stationary set, anything may happen with the entangle- 
ment for short times: it may die, resurrect, oscillate, etc. 
It is also important to mention that such analysis does 
not depend on the specific entanglement quantifier used 
to follow the dynamics, only the assumption that it is 
continuous and strictly positive on entangled states. 

Given a dynamics that fits in case ii) one can in gen- 



eral find examples of initial states whose entanglement 
die asymptotically or suddenly [ill]. An interesting way 
to have a global view of the properties of this dynam- 
ics on this respect is through the question: if one pick a 
random initial state, what is the most probable situation, 
asymptotic or sudden death? That is, if the dynamics can 
exhibit both of these properties, what is the most typi- 
cal one? To answer this question one must formulate it 
properly. Fixed a dynamics for a composite system with 
state space D with a suitable probability measure P on 
it and a continuous entanglement quantifier e : D M_|_ , 
with e{E) C (0, oo), we define the following events (sub- 
sets of D, in the language of probability theory) whose 
probabilities may be of interest: 

• States that exhibit sudden death of entanglement: 
SDE ^ {p e VN\3to,ti such that E{p{to)) > 
and E{p{t)) = for alH > ii}; 

• States that exhibit asymptotic death of entan- 
glement: ADE = {p £ VN\3{tn)?^=i,t„ 

oo, such that E{p{tn)) > and lim(^oo-E'(p(^)) = 

0}; 

where p{t) denotes the time t evolution of initial state p 
according to the dynamics. Note that these definitions do 
not coincide strictly with the common sense of such no- 
tions since in general one only looks for initial states that 
already have some entanglement, which is not necessary 
here: an initial separable state can, in principle, acquire 
some entanglement that will subsequently die (suddenly 
or asymptotically). The strict notion would be given by 
the events: 

• SDE' = SDE n E; 

• ADE' = ADE n E. 

If the dynamics exhibit asymptotic entangled states, one 
can also look to the events: 

• The states exhibit entanglement asymptotically: 

AE = {p e VN\3ta,c > where E{p{t)) > 
c for alH > <o}; 

• An initially separable state acquire entanglement 
asymptotically (sudden birth of entanglement): 

SEE = {p e Vn\E{p) = Oand3to,c > 
where E{p{t)) > c for all t > to}; 

(note that SEE ^ AEr\S). 

Instead of choosing a specific probability measure to 
deal with, our results will only require that it is non- 
singular, i.e., sets contained in sub-manifolds of D with 
dimensions strictly smaller than the dimension of D have 
zero probability. The problem of computing the proba- 
bility (or volume) of the event (set) 5" exactly is still 
an open issue for the most natural probability measures. 
Though, several bounds and estimates exist for several 
probability measures and events [l^ . [l7| . 
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III. EXPLICIT EXAMPLES 

Given the general picture we may look now to some 
concrete examples where most of them, as we will see, 
are very natural and experimentally feasible. The sim- 
plest type of dynamics, namely one that is convex-linear, 
Markovian and completely positive, will suffice to provide 
rich examples. It will be sufficient to work with the sim- 
plest composite system, two qubits, in order to exhibit 
examples for all classes of dynamics. Considering that we 
are dealing with a system with finite Hilbert space, we 
may apply Lindblad theorem and describe the map by 
an ordinary, linear, first order differential equation with 
the form given by [18'|: 

^ = C[p] = -^[H,p]+V[p], (la) 

where 

m = E - A^A.p - pA]A,), (lb) 

j 

7j are real constant numbers, Aj general linear operators, 
H a Hermitian operator. This type of dynamics have the 
advantage that it is simple to find its asymptotic states: 
in general one just have to look to the kernel of the "su- 
peroperator" C, which is a Hnear operator that can be 
understood to be defined over the set of 4 x 4 complex 
matrices or the subset of Hermitian matrices, a real vec- 
tor space. Of course, since the set will be given by the 
kernel of linear map, it is always given by the intersec- 
tion of a subspace of Hermitian matrices (the kernel of 
C) with the set of mixed states. It is courious to note 
that to find the two missed examples in Ref. [l3| we had 
to allow for non- autonomous Lindblad equations, that 
is, equations with the same form but with parameters 
7j varying in time. For this type of dynamics the set of 
stationary states do not need to be the intersection of a 
subspace with the set of quantum states. 



D (if it is grater than or equal to zero, respectively). The 
determinant of the partial transpose, on the other hand, 
gives us complete information about its entanglement be- 
cause it is known fl9|| that the state is entangled iff the 
determinant is strictly negative. Thus, this determinant 
tells us if the state is in the interior of the set of separable 
states (if it is grater than zero), in the border dS (if it is 
equal to zero) or inside the set of the entangled states (if 
it is strictly negative). At last, if p is a given state with 
d — Detp and (f — Detp'", we can have: 

i) c? > and cF > 0: the state is in the interior 
of D and in the interior of S relative to D, i.e., 
belongs to S — dS {e.g., the completely mixed state, 

Pmix — -^/4) 5 

ii) c? > and dF — 0: the state is in the interior of D 
and in dS (e.g., the state ^Pmix + kpaingiet, where 
Psingiet rcfcrs to the state in Eq. 0|) with a = 0, 
b^-c = 1/2); 

iii) d > and d^ < 0: the state is in the interior of 
D and belongs to E {e.g., the Werner states (20| 

PPmix 

+ (1 -p)psingiet, for < p < 2/3); 

iv) c? = and d^ > Q: the state is in the border of D 
and m S — dS (recording that we defined dS as the 
boundary relative to D, while dD is relative to A). 
For instance, if a > 6 > 0, 2a + 26 = 1, and |c| — b: 

/ a \ 
6 c 

^ = c 5 ' (2) 
\ a/ 

v) d — and d^ — the state is in dD n dS {e.g., a 
separable pure state); 

vi) c? = and d^ < the state is in the border of D 
and belongs to E {e.g., Psingiet)- 

With all these tools in hand, we can go to the examples. 



Localizing a two qubit state in the set of all states 

The set of all quantum states for a composite system 
can be divided geometrically according to the dichotomy 
{IntD, dD} and the trichotomy {IntS*, dS, E}. Dealing 
with the special case of two qubits has the advantage 
that one can easily infer the location of a state accord- 
ing to this subdivision with the help of Detp and Detp'", 
the determinants of the state and of its partial trans- 
pose. Both of these functions are continuous in all natu- 
ral metrics, Hilbert-Schmidt, etc, i.e, we know that small 
perturbations of a state in a given metric implies small 
perturbations of the values of both quantities. So if, e.g., 
both of them are positive for a given state, we can find 
a neighborhood of that state where these quantities re- 
main with the same sign. Then, the determinant of the 
operator tell us if it is in the interior or in the border of 



Case la): One asymptotic state in Int(S') 

Perhaps the most natural example of this situation is 
the case where both qubits, which we shall call A and B, 
are spatially well separated two level atoms, interacting 
with thermal fields. The separation between them im- 
plies that the thermal reservoirs are independent. The 
Lindblad equation that describes this dynamics is given 
by: 

^ = ^[Ha + Hb,p\ +Va® I[p\ +1® Vb[p\, (3a) 
where 

T^Ap] = 7i(2o-+,iPO--^i - a-^ia+.ip - pa-^ia+.i) + 

7-(2cr_,i/9CT+,i - a+^ia-^ip - pa+^ta-^i), (3b) 
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with cr±i being the Pauli operators for qubit i, 
Hi = —^az,i the Hamiltonian for qubit i and 7^,7^' are 

non-negative constants (related to the average photon 
number in the field, the atoms polarization, their cou- 
pling to the environment, etc.). 

It is easy to show that the system will evolve to a 
product state with both qubits in their respective Gibbs 
states, Z~^e~^^\ Zi = TTe~^^\ If the temperature 
is positive, the resulting state is a product state with a 
diagonal density matrix (in the product basis) with every 
diagonal entrance being non-zero. We then have that 
Pst = p^st ^nd that Detp'" = Detp > 0. As mentioned 
earlier, if an initial state have some entanglement it will 
certainly die at finite time. 

The events defined in sec. [ll] are trivial in this case: 
ADE = ADE' = SEE ^ AE = %, and SDE = SDE' = 
E, so P{SDE) = P{E) or, P{SDE\E) = 1, that is, the 
only condition for having entanglement sudden death is 
that the initial state is entangled. To calculate the exact 
probability of this event is thus as difficult as determining 
the volume of the set of separable states pL6l |. 



Case lb): Several asymptotic states in IntS 

To obtain an equation of motion for the state satis- 
fying this propriety, namely, being a relaxing dynamics 
with more than one asymptotic state but all of them in 
the interior of S, we had to appeal to a non-autonomous 
Lindblad equation. A dynamics that achieve the de- 
sired result would be given by a Lindblad equation with 
the same form as the one used in the last section, de- 
scribing two qubits interacting with independent reser- 
voirs, but now, with the coupHng "constants" decaying 
exponentially. That is, performing the correspondence 
7i ^ 7ioexp(— Ki). The physical situation correspond- 
ing to the equation, although artificial, is certainly not 
prohibited: in principle, one can have a good control of 
the interaction of the qubits with their reservoir and turn 
it off exponentially. 

To prove the result, let us write the dynamical equation 
in the form (in the interaction picture): 

^^p{t) = e--'V[p{t)], (4) 

where V is the dissipator of the Lindbladian in the last 
example. For p{t) a solution to this equation, we can 
define p[t) — {pog)(t), where 

g{t)= f e-^''dt' 
Jo 

is an invertible function. Substituting p in Eq. ^ we 
obtain an equation of motion for it: 

Im^vipit)]- (5) 



That is, p obeys the same dynamics of two qubits in 
independent thermal reservoirs with constant coupHng 
in time, with known solution. To find the asymptotic set 
for the dynamics of Eq. Q is sufficient to note that, since 
P{t) = iP°9~^)(t), then p{t ^ 00) ^ p{9~^{t 00)) = 
p{l/^). 

Geometrically, the autonomous dynamics given by 
Eq. 1(5]) deforms continually the set of states D to the 
point pGibbs (i.e., provides an homotopy between them), 
while the time varying version reparametrizes this defor- 
mation. The set of asymptotic states of Eq. ([4]) is then 
given by this deformation in the intermediate time k~^. 
Making k small enough, we can assure that the asymp- 
totic set is entirely contained in IntS', since pcibbs belongs 
to IntiS*, an open set. 

Of course, the events SDE, ADE, etc., and their re- 
spective probabilities, are exactly the same as in the last 
example. 

We note finally that, although the discussion about 
entanglement does not depend if we are dealing with 
the interaction or Schrodinger pictures (because the cor- 
respondence between them is given by local unitary 
transformations), the dynamics is not relaxing in the 
former. Since the state will be given by ps{t) = 
exp{iHt)p(t) exp{—iHt) and limj^oo p{t) will not, in gen- 
eral, commute with the exponentials, the state evolution 
Ps{t) will not converge. Nevertheless, the dynamics will 
have an asymptotic set in the general sense discussed in 
Sec. ini namely, although an initial state does not nec- 
essarily converges, one can find open sets such that the 
state trajectory will be confined inside them after a cer- 
tain instant of time. In this particular example, one can 
find such open sets that are entirely contained in S. 

Case 2a): One asymptotic state in dS 

Eqs. ([3]) also provides an example where we have only 
one stationary state in the border between separable and 
entangled states, namely, the case where the qubits are 
subjected to two independent thermal reservoirs at null 
temperature. In this case the stationary state is the pure 
state Pst — |00) (00|. Again, it is diagonal in the compu- 
tational basis so Detp^j = Detpst = 0. Then, a neigh- 
borhood of this state always contains separable as well 
as entangled states. As mentioned in Sec. m in this ex- 
ample, depending on the initial state, both behaviors can 
happen: asymptotic and sudden death of entanglement. 
In fact, given an initial state with matrix elements pij one 
can shown that the determinant of the partial transpose 
of the state in time t will be given by: 

Detp^ (t) = e-'^'^'Det [p' + p" {t)\ , (6a) 

where 

(Pii P\i P\-i P23 \ 

^"t^^^ ^*'lo^i' , (6b) 

Pl3 Pl4 Pll + P33 P34 + 2Pi2 i 

P23 P24 + 2pt3 P34 + 2pi2 1 / 
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and p"{t) is a matrix which depends on p but where all 
elements decay (exponentially) to zero. Hence, as long 
as Det/c' ^ 0, the asymptotic sign of Detp'"(t) will be 
given by the sign of the determinant of p' . By assuming 
non-singular probability measure in the set of quantum 
states, we conclude that the event defined by the con- 
dition Detp' = has zero probability and can be dis- 
carded to compute the probabilities of ADE{= ADE') 
or SDE{= SDE'). From the form of p' it is easy 
to find initial states such that Detp' is strictly less or 
strictly greater than zero, so small balls (with positive 
probability) around these states also have the same sign 
for this determinant. As a consequence, we have that 
P{SDE) > 0, P{ADE) > 0, the actual values depend 
on the specific measure used. The point is, with no 
additional requirement on the measure, both situations, 
asymptotic or sudden death, can be found for this dy- 
namics. 

Since this dynamics do not have asymptotic entangled 
states one have SEE = AE = %. 



Case 2b): More than one asymptotic states with 
points in the border of S with E 

For more than one asymptotic state in this geometric 
situation we can distinguish four subcases, as discussed 
below. 

All other points belong to \ntS . Two non-interacting 
qubits subjected to two independent phase reservoirs pro- 
vide an example. The dynamics (interaction picture im- 
plied) is given by: 



dp 
lit 



Va^I[p]+I(S>Vb[p] 



where 



T^dp] = l{o'zAPO-z,i - p), 



(7a) 



(7b) 



with 7 a positive constant. This dynamics may be im- 
plemented experimentally for ions in a trap 1.11] . The 
reservoir would be given by applying z-directed magnetic 
fields with random and independent magnitudes on each 
ion (m (the qubits encoded in the electronic spin of the 
ions) . It is easy to show that if we write the initial state 
in the computational basis the evolution will be given by 
exponential decays of all non-diagonal terms and all the 
diagonal ones will remain constant. So the set of asymp- 
totic states will be given by the three real parameters set 
(an intersection of a four dimensional subspace of the set 
of Hermitian matrices with the set of states): 



Pst 



Pi 

P2 

P3 

p4 



(8) 



In this case we have that all asymptotic states are 
diagonal in the computational basis and again we have 
Det/o —Betp^. Two situations are possible: these deter- 
minants are zero or positive. Again, entanglement can 
die asymptotically or suddenly as the following initial 
states illustrate: 



Pit = 0) = 



/ Pi 

P2 c 

c p3 

V Pi 



(9) 



with |c| > (as a consequence, P2 > and P3 > 0). The 
evolution will be given by states with the same form but 
with \c{t)\ decaying exponentially, so dF{t) = P2P^[piPA — 
|c(i)p). Then, it is evident that, if pi or pA are initially 
zero, entanglement will decay only asymptotically to a 
state in the border of S. But if both of them are non- 
zero and PiPa < |c(0)p, then it will die suddenly while 
the state converges to (a state in) the interior of S. For 
PiPa > |c(0)P the complete trajectory will remain in S. 

Although examples of both situations can be pro- 
vided, the typical case is definitely sudden death of 
entanglement [20|- As this dynamics do not exhibit 
asymptotic states with entanglement, SEE = AE = % 

All other points belong to E. For this case we chose 
a situation where both qubits are identical (but distin- 
guishable) and interact collectively with a common reser- 
voir, as it happens with two spatially close two level 
atoms (close compared to the wavelength defined by their 
transition) in a thermal field. The dynamics of this situ- 
ation can be described by the following master equation 
(also in the interaction picture) [i^: 



dp 
'dt 



7(2J_pJ+ — J^J^p — pJ^J^) 
+-f'{2J+pJ_ - J_J+p - pJ-J+), (10a) 



with J± = <J±^A + <7±,B- A convenient way to analyze 
this dynamics is to write the equations of motion for the 
density matrix elements in the basis composed by the 
states {|11) , 1^'+) , |00) , |*-)}, resulting: 



pii = -27P11 + 27'p22, 

P22 = 27(pii - P22) + 27'(p33 - P22), 
P33 = 27P22 - 27'p33, 
P44 = 0, 

P12 = -27P12 + 27'p23 - 7'pi2, 

Pl3 = -7Pl3 - 7V13, 
Pl4 = -7Pl4, 

P23 = -7P23 + 27P12 - 27'p23, 



(10b) 



P24 = -7P24 - 7 P24, 



P34 



-7 P34. 



with Pi > for i = 1, 4 and J^i 



Pi 



1. 



The reservoir at zero temperature corresponds to the 
case 7' = 0. It is easy to see from the equations of motion 
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that the complete subspace span {|00) , Ivf")} is station- 
ary under this dynamics. By convexity, the stationary 
states have the following form: 



Pst 







1 - P44 P34 

P34 P44 



(11) 



and can be identified with a Bloch ball inside D. All 
states have null determinant, so all of them are at 
the boundary of -D. It is readily seen (representing 
these states in the computational basis) that Detp'" = 
— (^44/2)^, which is null only if P44 = and is negative 
otherwise, so the set do not have any points in the inte- 
rior of S: this Bloch ball just touches the set of separable 
states in one point. Some things can be inferred immedi- 
ately from the geometry of this set: a) the entanglement 
of the system may never die (the singlet state is station- 
ary, for instance); b) it can be created: take any initially 
separable state p with non-zero population in the singlet 
state; c) In principle, the entanglement can die asymp- 
totically or suddenly. In fact, initial states leading to this 
situation exists but only a) and b) are "typical". 

A helpful fact about this problem is that the singlet 
population is constant through the evolution so, if this 
population is positive on the initial state, it will con- 
verge to an entangled state. Since the event formed by 
all states with non-zero population have probability one, 
we immediately infer: P{AE) = 1,P{SBE) = P{S), or 
P{SBE\S) = 1, that is, if one chooses randomly an initial 
state, regardless if it is entangled or not, it will evolve to 
an entangled state with probability one. From this we im- 
mediately see that P{ADE) = P{SDE) = P{ADE') = 
P{SDE') = 0, i.e, the probability to choose an initially 
entangled state whose entanglement will vanish is zero. 

Nevertheless, one can find atypical specific examples 
exhibiting SDE and ADE. Consider, for instance, the 
family of initial states where the only non-vanishing 
matrix elements (in the basis mentioned above) are 
Pii, P22, P33- From Eqs. IjlObp it follows that those will 
continue to be the only non- vanishing elements. If 
also pii = their behavior is quite simple: piiit) = 

0,P22(t) = P22e-2^*,P33(t) = 1 - P22{t). So if p22 ^ 

the state will remain entangled for all times (mixture 
of a Bell state with an orthogonal separable state) and 
will die asymptotically, i.e., exhibit ADE. On the other 
hand, if pu ^ the behavior of these matrix ele- 
ments is still simple and the determinant of the partial 
transpose will acquire the following form: Detp^(t) = 
piie~'^'^* + P{t)e~^'^* , where P{t) is a second degree poly- 
nomial with coefficients determined by the initial den- 
sity matrix elements. Since pn ^ 0, this determinant 
will be positive after a certain instant of time, i.e., the 
state will be always separable after that instant. If, e.g., 
P33 — 0, Pii 7^ 0, P22 7^ the initial state is entangled 
and therefore will exhibit SDE. 

Some points belong to IntS and others to E. The 
reservoir used in the last subcase, if taken at positive 



temperature, provides this example and, to simplify the 
problem we take the infinite temperature limit (7 = 7' in 
Eq. (|10ap ). It is interesting that, irrespectively of temper- 
ature, the singlet state is stationary and also the singlet 
population of any state (the singlet spans a one dimen- 
sional decoherence free subspace for this model). From 
the equations of motion immediately follows that the sta- 
tionary states are: 



Pst 



V 



^ \ 



^ 

p J 



(12) 



where p is the singlet population of the state. That is, 
they are the Werner states (with a different parametriza- 
tion). 

The determinant of the partial transpose (with re- 
spect to the computational basis, of course) is simply 
(3 — 12p^)/36 being negative only if p > 1/2. The set 
of stationary states forms a line segment in D with both 
ends, those with p = or p = 1, on the border of D, one 
of them in the interior of S (relative to D) and the other 
in E, respectively, and the line intersecting the border 
between S and D when p = 1/2 (see Fig. [2]). 




jDmix triplet 



Psinglet 



FIG. 2: Set of asymptotic states for two qubits interact- 
ing with a common reservoir at infinite temperature. Here, 

P™x tripUt = Kill) (111 + |*+> {*+| + |00> (001). 

Since the singlet population remains fixed in the dy- 
namics it allows us to identify the asymptotic state of any 
given initial condition. An initial state will have non-zero 
entanglement asymptotically if, and only if, ^44 > 1/2 so 
we have P{AE) = P{D>i/2 = {p e D\p44 > 1/2}) > 0. 
Of course P{SBE) — P{D^i/2(^S). Since a state can ex- 
hibit ADE iff it relaxes to a state in the border between S 
and E, we have P[ADE') < P{ADE) < P{ADE Ci {p e 
D\p44 — 1/2}) — 0. So ADE is atypical for this dynam- 
ics but SDE, on the other hand, have a non-zero proba- 
bility. In fact, an initially entangled state have SDE iff 
P44 < 1/2, so P{SDE') = P{En{p e D|p44 < 1/2}) > 0. 

All points belong to dS. The combination of two 
reservoirs used in former examples will provide this case. 
If we have qubit A subjected to spontaneous decay and 
B to a phase reservoir the system will have the desired 
behavior, a situation that may occur experimentally if we 
entangle an atom in vacuum with a spin subjected to a 
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stochastic magnetic field. That is, the system dynamics 
would be described by a master equation of the form l(7a|) 
(again in the interaction picture), but with Va given by 
Eq. lUEl (with i = A and -j'a = 0) and by Eq. ^ 
(with i — B). It is easy to see that the set of asymptotic 
states will be constituted by the product states where A 
is in the |0) state and i? in a state described by a diagonal 
matrix (in the computational basis), so the global states 
reads: 



Pst 





p 



1 -p , 



(13) 



for < p < 1. Whatever the value of p we have 
Detpst =Detp^j = so they indeed belong to dS. Again 
we have SBE — AE — for this dynamics, since there 
are no entangled asymptotic states, but to analise the 
probability of the other events we use the exact solution 
for the dynamics and write the determinant of the partial 
transpose in the form: 



Detp(0^ = /i(p) 



fn{p)e 



-A„t 



(14) 



Pll 


Pl2 








Pl2 


Pll + P22 














P33 


P34* 








P34 


P33 + P44 



where the functions depend on the initial state only, 
while Ai < A2 < ... < A„. In this way, as long as /i 7^ 0, 
the asymptotic sign of Detp(t)'" will by given by the sign 
of /i. Denoting by 7^ and 7b the decay rate for each 
reservoir, it so happens that Ai = 2'~fA and /i = Detp', 
where: 



(15) 



Since this matrix is positive definite (given that p is), 
the system will reach its asymptotic state from the in- 
terior of the separables if fi{p') > 0. But the event 
fi{p') — have zero probability, so we may conclude 
that P{SDE) = P{SDE') = P{E), while P{ADE) = 
P{ADE') — 0, that is, a sorted initial entangled state 
will exhibit sudden death of entanglement with certainty, 
in contrast with case 2a) where sudden and asymptotic 
death both had positive probabilities. Still, it is possi- 
ble to find specific states where asymptotic death takes 
place. For instance, consider the set of initial states: 



/O 

P22 P2-i 

P23 P33 

\ 0, 



(16) 



Case 3a): One asymptotic state in E 

The most natural way to realize a dynamics with this 
property is through a thermal reservoir. This time, 
though, interacting qubits and a common thermal reser- 
voir are needed, that is, a reservoir that take any ini- 
tial state to the Gibbs state exp(— /?i7), with Z = 
Trexp(— /3i/) and H stands for the Hamiltonian describ- 
ing the closed dynamics of the qubits. Typically, the 
ground state of interacting qubits Hamiltonian is non- 
degenerate and entangled, so if /3 is large enough we ob- 
tain the desired dynamics. 

Dynamics with these asymptotic states can be engi- 
neered using Lindblad autonomous equations, at least 
formally. Actually, fixed an arbitrary state for the sys- 
tem, there are many Linbdbladians that have this state 
as the only asymptotic state, in particular there are ones 
with only one Lindblad operator and null Hamiltonian 
part jl^l. The specific Lindbladian of course, will depend 
on the specific interaction between qubits and reservoir 
(if the dynamics could be described by a Lindblad equa- 
tion in the first place). 

To give a more specific picture, consider, for instance, 
two interacting qubits described by the following Hamil- 
tonian: 

H = \<ja^,A + \<^(J^^B + 5(0-+, AO--, B + cr-,AO'+,B) (17) 

with g positive constants satisfying g > uj (i.e., strong 
coupling limit). The eigenvalues for this Hamiltonian 
are, in crescent order, — g, — w, uj, g, with respective eigen- 
vectors l^*-) , |00) , |11) , 1^'+), leading to an entangled 
ground state. Denote by \i), i = 1, . . . ,4 these eigenvec- 
tors according to their eigenvalues order. We may con- 
sider a thermal reservoir at null temperature that induces 
decays between any two of these states in a Markovian 
way, such that the dissipator would be: 

where (Xy = \i) (j| and 7^^ are non-negative constants. A 
dissipator of this type can be derived from a microscopic 
model, for instance, adapting the results of Ref. [23| to 
the Hamiltonian considered here . 

As in case la), the events and probabilities we are in- 
terested in are trivial: SBE — S,AE = D, i.e., every 
initial state will acquire entanglement for large times, 
in particular the separable ones, so P{SBE) = P{S) 
and P{AE) = 1. Since the entanglement never vanishes, 
ADE = ADE' = SDE = SDE' = 



Case 3b): Several asymptotic states in E 



The partial transpose determinant will be then 
T>etp{t)^ = -|p23(i)PP22(i)P33(i), being negative for all 
t if P23,P22 and P33 are initially different from zero, so 
the entanglement dies asymptotically. 



Examples for this case can be provided just by the same 
trick used in case lb): we take any Lindbladian with only 
one asymptotic entangled state and insert a time vary- 
ing coupling which multiplies the dissipator. The same 
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reasoning can be applied with respect to the asymptotic 
states for the subsequent dynamics (in the interaction 
picture), so, if the decay rate of the couphng is small 
enough, the set of asymptotic states will be constituted 
by a small "blurring" around the asymptotic state of the 
dynamics with constant coupling. 

Contrary to case lb), though, in what entanglement is 
concerned, it is important now whether the dynamics is 
given in the Schrodinger or interaction pictures, because 
their correspondence is given by global unitary transfor- 
mations. By the same reason as before, the dynamics will 
not be relaxing in the Schrodinger picture, but one can 
still find a non-trivial asymptotic set, this time, entirely 
contained in E. Indeed, diminishing the decay rate of the 
reservoir coupHngs, we can diminish at will the diameter 
of the set of stationary states in the interaction picture 
which, by its turn, always contain the Gibbs state of the 
system. Now, unitary transformations are isometrics for 
practically all relevant metrics, so the set of asymptotic 
states in the interaction picture is mapped to sets with 
the same diameter in the Schrodinger picture. But these 
unitary transformations have the Gibbs state as a fixed 
point, hence these sets always contains it. Since E is 
open, given that their diameter is small enough, we can 
be sure that they always fall entirely inside of it. 

As a consequence of the discussion in the above para- 
graph, the events and probabilities we are considering in 
this paper are identical to the ones in the last example. 



IV. CONCLUSIONS 

In this paper we review the classification of the possible 
dynamics of entanglement based on the relative geome- 
try of the sets of asymptotic and separable states. We 
provided examples for all possible classes, including the 
previously unknown cases with more than one asymp- 
totic state, but avoiding the boundary dS. In giving 
those examples it was sufficient to use two-qubit dynam- 
ics dictated by equations of motion in Lindblad form 
(including non-autonomous dynamics exactly for those 
previously hard examples). In each case, the existence 
of sudden death of entanglement, asymptotic death of 
entanglement, sudden birth of entanglement and asymp- 
totic entanglement were analyzed from a more precise 
point of view, looking at the probabilities that each of 
these phenomena occur if one choose a random initial 
state and a suitable probability measure on the set of 
quantum states. 
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